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Abstract 

in , 

04 ' We consider the movement of a particle advected by a random flow of the form v + (5F(x), 

with V € M.'^ a constant drift, F(x) - the fluctuation - given by a zero mean, stationary random 

Q field and 5 ^ 1 so that the drift dominates over the fluctuation. The two-point correlation matrix 
■r>^„\ „f 4-1 I c„lJ J l-.l2a-2 I — I , -^x. „ ^ i rrii _ t^.i „ r i_ <•„_ 



R(x) of the random field decays as |xp"~^, as |x| +oo with a < 1. The Kubo formula for 
the effective diffusion coefficient obtained in |17| for rapidly decorrelating fields diverges when 
^1 1/2 < a < 1. We show formally that on the time scale the deviation of the trajectory 

from its mean y{t) = x(i) — -vt converges to a fractional Brownian motion Ba{t) in this range 
of the exponent a. We also prove rigorously upper and lower bounds which show that E[|y(i)p] 
converges to zero for times t <C (5^^^" and to infinity on time scales t ^ as (5 ^ when 

a e (1/2, 1). On the other hand, when a < 1/2 non-trivial behavior is observed on the time-scale 

o 
i> 
o 

\ The position of a tracer in a random flow is described by the ordinary differential equation 



1 Introduction 



^5|^=V(t,X(t;x)), X(0;x)=x, (1.1) 

. at 

^ I where V(t, x) is a random field. This random model is frequently used to describe the motion of 

a particle in a turbulent flow of fluid, where V(t,x) is the Eulerian velocity field. The long-time 
behavior of solutions of ()1.1() has been extensively studied, especially when V(i, x) is sufficiently 
\ rapidly mixing in time and space ~ see j21j for an extensive overview. Roughly speaking, one 

expects that the trajectory X(t; x) behaves diffusively if the velocity field decorrelates sufficiently 
fast in time and space while an anomalous behavior is observed if V(t, x) has long range correlations. 
This problem remains largely open when V(x) is both time-independent and of zero mean, as then 
there is no mechanism to "move the particle around" and employ the decorrelating properties of the 
random field. Therefore, to simplify matters, we assume that the flow satisfies the Taylor hypothesis, 
that is, it has a non-zero mean drift which dominates over the amplitude of its fluctuations. This 
means that we can write V(x) = v -|- (5F(x), where v 7^ is a constant vector with, say, |v| = 1, 
and the parameter (5 <C 1. The fluctuation F(x) is spatially homogeneous, of zero mean and is 
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divergence-free: V • F(x) = 0. It has been shown in ^Jj that if the field F(x) is sufficiently strongly 
mixing in space then the scaled trajectories 



y(t/(52;x) :=x(t/52;x) - vt/(52, t>0 (1.2) 

converge, as (5 | 0, as continuous stochastic processes, to a zero mean Brownian motion starting at 
X with the covariance matrix D = [Dij] given by the Kubo- Taylor formula 

Dij = ^ j [Rijivt)+Rjiivt)]dt, i,j =,..., d. (1.3) 



Here R(x) = [i?jj(x)] is the covariance matrix of the field F(x), that is, Rij{x.) = E[Fi{x.)Fj{0)]. 

In the present paper we are interested in the situation when the correlation tensor decays slowly 
in space so that the diffusion matrix given by the Kubo- Taylor formula is infinite and the standard 
diffusion limit may fail. More precisely, we consider Gaussian fields which are isotropic and satisfy 
locally the self-similarity hypothesis. This means that the two-point correlation tensor satisfies 

Rij{^)= [ e'''-^R,,{k)dk, (1.4) 

with the power spectrum 

^^^•(k) = ^i^!^r,,(k), i,j =,..., d. (1.5) 

Here a(|k|) > is a non-negative bounded, measurable function, supported in a finite ball {|k| < K} 
for some K > 0, continuous at k = and with a(0) = 1. The factor Tij(k) := 5ij — kikj, where 
k = (ki, . . . , kd) ■= k/|k|, ensures incompressibility of the flow. To guarantee the integrability of the 
spectrum we assume that a < 1. The rate of decay of the correlations of the field is then given by 

R^j{^) ~ |x|2°~^ for |x| » 1. (1.6) 

A simple calculation shows that in this situation the diffusion matrix given by the Kubo formula 
(|1.3j) is infinite, provided that 1/2 < a < 1. The goal of this paper is to find the proper time scaling 
for the tracer trajectory when the Kubo- Taylor formula diverges. In other words we are looking for 
H such that for any p > 

^lim^E|y(t/52^(^+'');x)|2 = +0O and lim_E\y{t/6'^^^^-P^ = 0. (1.7) 

The result of Kesten and Papanicolaou, see Theorem 4 of ^7], guarantees that for R(x), which is 
decaying sufficiently fast, or, equivalently, for a < —N with a large enough > one has H = 1. 
Note, however, that the Kubo-Taylor formula ()1.3|) itself makes sense for any a < 1/2, that is, in 
a much larger range of a. We shall prove rigorously in Section El below that then indeed H = 1 for 
a < 1/2, with a slightly modified definition of the limit appearing in condition (|1.7|) . 

According to ()1.6() . increasing the value of a leads to strengthening of the correlations of the 
corresponding field (and degraded mixing properties). We shall prove, see Theorems 13.11 and 13.21 
below, that the level a = 1/2 is critical in the sense that the value of the exponent H changes at 
a = 1/2 from H = 1/2 for q < 1/2 to = l/(2a) < 1 for 1/2 < a < 1. 

We also present a formal argument which provides more refined information on the process y{t) 
for a G (1/2, 1). In addition to identifying the time scale t ~ 0{6~'^^), H = 1/(2q) on which y{t) 
behaves non-trivially, we are able to compute formally the limit of the process 

ysit) = y(t/<52^;x) = x(V5'^;x) - vV5'^. 
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as ^ 0. It turns out that the hmit of ys{t) is a super-diffusive fractional Brownian motion Ba{t). 
This should be contrasted with a standard Brownian motion limit for the process y(t/5^) obtained in 
jl7j for rapidly decorrelating fields F(x). We recall that a Gaussian, continuous trajectory, R'^-valued 
process {Ba(t))t>o is called a fractional Brownian motion (FBM) with Hurst exponent a G (0, 1) and 
the covariance matrix D if for any t,h >0 the (normal) random vector Ba{t + h) — Ba{t) has mean 
zero and the covariance D/i2°. In particular, E(S„(t) (g) = Dt^" and the tracer propagation 

is close to ballistic as a | 1 and correlations persist at infinity. Note that for a = 1/2 the FBM 
becomes the "usual" Brownian motion matching the result of ^7] for rapidly decayng correlations. 

Anomalous diffusive behavior of passive scalar s has been studied extensively in the physics lit- 
erature - we mention [HI E3 1213 ESI I2E] without any attempt at completeness. However, the 
majority of physical papers consider time-dependent (typically white in time) random flows with a 
non-zero diffusivity. Super-diffusive limits for particles in a random flow have been also obtained 
rigorously in various situations ~ for instance, in j^lEllllISI for the random shear layer flows, in 
for time-dependent random flows, in ^Hl for a mean zero, random flow with a positive molecular 
diffusivity and in jTlElj for multiscale periodic flows. The novel aspect of the present paper is that 
the flow is time-independent and the molecular diffusivity is equal to zero. 

The paper is organized as follows. The aforementioned formal derivation of the fractional Brow- 
nian motion limit for a G (1/2,1) is presented in Section [21 We describe in Section 01 upper and 
lower bounds for the mean square displacement of trajectories that show the super-diffusive scaling 
for Q G (1/2,1) and diffusive behavior for a < 1/2. The main results are Theorems 13.11 and 13.21 - 
they are proved in Sections and IHI Some auxiliary bounds are proved in Section IHl Appendix 1X1 
contains a brief review of multiple stochastic integration. 

Acknowledgment. LR was supported in part by an Alfred P. Sloan fellowship. This work has 
been also supported by ONR and NSF grant DMS-0604687. 



2 A formal limit for the advection equation 
The advection equation 

We present here a formal computation that leads to a FBM limit for the solution of the advection 
equation 

^ + (v + 5F(x)) • Vx05 = 0, (/<5(0,x) =no(x). (2.1) 

Here v G M'^ is a fixed mean flow with |v| = 1 and F(x) is a spatially homogeneous Gaussian random 
field with the two-point correlation tensor [i?jj(x)] of the form (|1.4|) and the power spectrum as in 
H1.5() . We can write then F(x) = f e*^ ''F(dk), where the spectral measure /^(dk) satisfies 



E 



Fn{dk)F^{dp) = Rnm{k)6{k + p)dkdp. (2.2) 



We will assume in this section that the exponent a G (1/2, 1) as this is the range where we expect 
the FBM behavior for a rescaled process in the limit S ^ 0. Let us set H = l/(2a) and introduce 
the rescaled time t' = and the moving frame x' = x -|- vt as well as the new unknown function 
us{t',x') = (j)s ^pTj-ix' — VpTj-^ Then ()2.1() becomes in the new variables 

dt' 

For convenience we introduce G(t,x) = — F(x — vt) and drop the primes, arriving at 



^ + 5i-2Hp|^^/_^_i_^ .^^,^^^0^ n5(0,x') = no(xO. (2.3) 



^ - S'-'^'G (^,x) • V^us = 0, us{0,x) = no(x). 



(2.4) 
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The main goal of this section is to argue formahy that E [n5(t, x)] converges, as (5 — > 0, to 

n(t,x) =E[Mo(x + 5D(t))]. (2.5) 
Here Bo{t) is the FBM with exponent a and covariance matrix D = [i^ij], where 



1 f e 



A. = ^ y ^^OT2r..(k)dk. (2.6) 

We will establish this limit by evaluating term-wise the expectation of the formal infinite Duhamel 
expansion for the solution of (|2.4|) and computing the limit of the main term. We will neither attempt 
to justify the expansion, nor try to estimate the error produced by the terms in the expansion which 
are formally of a smaller order. The main difference with the corresponding calculation in the 
case when correlations are rapidly decaying (the Brownian motion limit case) is that all Feynman 
diagrams contribute to the limit and not only the so-called ladder diagrams. 

An infinite expansion for the solution 
We rewrite (|2.4|1 as an integral equation: 

t 

U5{t,x) = uo{x)+5^-^^ j G (^,x) • Vxn5(si,x)(isi. (2.7) 



Iterating 1)2. 7(1 by substituting for ug{si, •) appearing on the right side we obtain a formal expansion 

+ 00 „ 

^,(t,x) = j; 5(1-2^)" / GH(^,...,^,x)dsi...d.. (2.8) 

Here, we have set G*-''-' (x) := no(x) and, assuming that G*-"^ (si, . . . , Sn,x) has been already defined 
for some n > 1, we let 

G("+^) (Sl, . . . , Sn, Sn+l, x) := G(s„+i, x) • VxG(") (Si, . . . , Sn, x) . 

The time integration region appearing in ()2.8|) is the simplex 

A„(t) := [{si, ...,Sn)eW :0<sn<...<si<t]. 

Passing to the Fourier transform in the definition of G^"^ and using induction we arrive at an explicit 
expression 

(^^] •F(dk, 

p=i 



M^o)dko. (2.9) 



G(")(si, . . . , s„,x) = /"e^iv.{kisi+...+k„s„)gix.(ko+...+k„) -Q 

P=i 

Next, we take formally the expectation of the infinite series in (|2.8j) term-wise and obtain 

+ 00 

E[us{t,x)]=J2U^){t,^), (2-10) 

n=0 

where 

„ 2n 

/„(5)(t,x)= 5^(^-2^)" J EG(2")(^,...,^,x)nd^.- 

We have used above the fact that the expectation of a product of an odd number of centered Gaussian 
variables equals zero. The next step is to evaluate each of the terms /n((5). 
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The main term in the expansion 

We will compute the individual terms 1^(5) using Feynman diagrams. Recall that given Gaussian 
random variables li, . . . , l2n the expectation E (Y1I2 • • • ^2n) is 



E ...Y2n)= n ^^(^P^^) 



(2.11) 



Here '^{n) is the set of all permutations (complete Feynman diagrams) of the elements of {1, ... , 2n} 
such that T o T = Id, and J-{k) 7^ k for all 1 < < 2n. The notation pij £ means that J-{p) = q 
(or that p and q are connected by an edge in the Feynman diagram J^) and p < q. 
Using (|2.11|) in the expression for /n(5) we obtain from (|2.9j) and (|2.in|) that 

In{6){t, X) = (_l)«52(l-2H)n ^ ^ J J ^-^.ik,s,S-^« +...+k,^s,^5-^^) 



A2„(i) 



9-1 



2n 



After evaluating the expectation above, using (|2.2|) . this expression becomes 



/„(<5)(t,x)=(-i)«<^2(i-2^)" 5] /•••/ n {«-p{--kp^} 

p— 1 q—1 2n 

}e^'^«-"no(ko)dkon^^- 



(2.12) 



The leading order term in (|2.12|) is 



/°(5)(t,x) = (-i)"52(i-2H)n ^ E /• • •/ n K • kp^l 

ii,...,i2n=l.FG5{n) A2„{t) P*?^-^ 



2n 



'~ko,ipko,i^Ripi^ (kp) dkp e*'^'^ '"'n (ko ) dko J| dsi 



(2.13) 



(_l)n^2(l-2H)„ ^ /•••/ n 



R 



i=l 

Sn S 



^27r^v ) ko • ko 



2n 



^ko)(iko JJdsj. 



2=1 



The reason why /^((5)(t, x) is indeed the main contribution to /„(5)(t, x) is as follows. Note that the 
two-point correlation function i?jj(x) decays algebraically for large |x|: 



i2ij(vs) = Aj(s)s'"~', 
with the matrix Dij{s), which converges, as s — > +00, to 



1 



rij{k)e'^-''dk 



l]^l2a+d-2 



(2.14) 



(2.15) 
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(a more explicit expression for Dij using Legendre polynomials may be obtained using the Hecke- 
Funk theorem (see e.g. jl5j . p. 181) but we will not need it here). Using the above information in 
(|2.i:-{|) we obtain: 



S2H 



2a-2 2n 

Y\dsi + o{l). 

1=1 



where D = [Dij]. However, and this is the crucial point in this calculation, as H = l/(2a) we see 
that 

2(1 - 2H)n - 2Hn{2a - 2) = 
and the powers of 5 exactly cancel each other. We conclude that /^(5)(t,x) = Jn{t,'x.) + o(l), with 

/„ 2n 
.../ e*^«-^(Dko • ko)" n \sg-Spf''-^uo{^o)dkol{ds, 



i=l 



In particular, /^(5)(t,x) is of order 0(1), as 5 ^ 0. On the other hand, all the other terms in 
(|2.12|) lead to expressions similar to (|2.1cil) but with Rij{x) replaced by its (possibly higher order) 
derivatives with respect to some of the spatial variables. However, derivatives of i?jj(x) decay faster 
than |xp~^° for large |x| - hence these terms produce a too high power of (5 as a factor, and vanish 
(at least, term-wise) in the limit (5 — > 0. 



Interpretation in terms of the fractional Brownian motion 

It remains now to relate Jn(t,x) to the fractional Brownian motion and sum all these terms. Note 
that the function 

/(si, . . . ,S2n) := ^ n I'^P ~ 

is symmetric in all of its arguments, that is, /(si, . . . , S2n) = /(^^(i)) • • • , s^(2n.))) where vr is an 
arbitrary permutation of {1, 2. . . . , n}. Using this fact we can rewrite Jn(i, x) in the form 



J„(t,x) 



E 



(2n)! — JO JO 



n 



l2o-2 



2n „ 

J{ds, j e^''«-"(Dko,ko)"no(ko)dko. (2.16) 



A simple but useful observation is that 



\2a-2 _ -2x1 



2^ /»CXD Q^kpSp 



n 



jj^w{dkp) 



where w{dki), . . . ,w{dk2n) are independent Gaussian white noises and Cq, > is given by 

/r(2a - l)sin(7ra)^ 



Cfy .- 



vr 



This follows from the fact that 



E 



|A;i|°-i/2 



w{dki 



oo \K2 



ko\'^-y 



^YHw{dk2) 



oo gifci (s—r) 
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for 



\2a-l 



dk = 2 



COS k 



dk 



TT 



fcSa-i r(2a - 1) sin(7rQ) ' 



The last equality follows e.g. from 3), [539] of p^. Hence, J„(t,x) has a representation 



J„(t,x) 



-1) 



n ft 



{2n)\ 



E 



2n 

n 

p=i 



—[J^w{dkp) 



2n 



Wdsi 

i=l 



iko-xi \2n. 



Dko,ko)"no(ko)(iko. 



Performing now the integrations with respect to Sj on the right side we conclude that 

2n 



(-1)' 

J„(t,x) = 



(2n)! 



n 

p=i 



c,P"(Dko,ko)"no(ko)e*'^°'^dko. 



iknx. 



Next, using the harmonizable representation of the standard fractional Brownian motion, see Propo- 
sition 7.2.8, p. 328 of [2lj, we obtain 



Jn{t,x) 



Here Ba{t) is a fractional Brownian motion with the Hurst exponent a and 

1/2 / ^ \ 1/2 



c : = 



TT 



TT 



Q!r(2Q) sin(a7r) J \2a^r(2a — 1) sin(a7r) 

Setting h := cca = l/(a\/2) we may now re-write J„ as 



(-l)"(2n- 1) 



n! 



(2n)! 

ikox 



2+2« 



-(Dko,ko) 



uo(ko)(iko. 



Coming back to (|2.1U|) we see that, as 5 ^ 0, 



E?i5(t,x) ^ ^(tjx) = ^ Jn{t,y^) = J exp jiko • x ^(Dko,ko) | ■uo(ko)<iko 

= y exp{i(x + BD(t))-ko}uo(ko)dko =E[no(x + BD(t))], 

where Bd (t) is a d-dimensional fractional Brownian motion with the exponent a and the covariance 
matrix D = D/(2a^). Therefore, we have formally established that H2.5() gives the limit of En5(t,x). 
In terms of the characteristics 



dX(t;x) 



v + 5F(X(t;x)), X(0;x)=x, 



(2.17) 

for the original advection problem (|2.1j) we have formally argued that the one dimensional statistics 
of the process 



t 



y(t;x) = X( ^;x) - 



g2H' 



with H = l/(2a), converge, as (5 — > 0, to the corresponding statistics of a fractional Brownian 
motion with the exponent a and diffusion matrix D given by (|2.6j) . In the next section we will 
show rigorously that indeed the process y(t) behaves trivially (either vanishes or E[|y(i)p] tends to 
infinity) on all time scales but t ~ 0(6^'^^). 
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3 The main results: rigorous bounds for the trajectories 



Preliminaries 

The probability space. Suppose that m is a positive integer and 'i?p(x) := (1 + |xp)~'', x G M*^, 
where p > d/2. Let Q be the Hilbert space of d-dimensional incompressible vector fields that is the 
completion of the space C^^.^ := {io G C^{R'^; R"^) :V^-uj = 0} with respect to the norm 

||a;||^ := J (|a;(x)|2 + \VM^)\' + ■■■ + | VXx)|2),?,(x) dx. 

We shall assume that m > d/2 + 1 so that any u; G $7 is of class of regularity by the Sobolev 
embedding theorem. 

The random field. The random field is set to be simply F(x;u;) := cj(x). Denote also F(a;) = 
{Fi{uj), . . . ,Fd{uj)) := uj{0). We suppose that P is a Borel measure given on Q, that satisfies the 
following hypotheses: 

(HI) it is Gaussian, that is, for any > 1, xi,...,XAr G M'^ we have F(xi), . . . , F(xjv) is a 
Gaussian, A^d-dimensional random vector, 

(H2) it is centered and homogeneous, that is, F is of mean zero and the two-point correlation 
matrix depends only on the relative position of the points: 

j Fi{z;uj)Fj{y;oj)dFiuj) = Rij{z - y), 

where R(x) = [Rij(x.)] is given by H1.4|) . As the measure P is Gaussian, this condition guarantees 
that P is invariant with respect to any spatial shift transformation Tx : f2 ^ $7, x G M'^ defined 
by Tx'^(z) := w(x + z). The existence of such a measure on Q is guaranteed e.g. by the results of 
|12j . see Section 2.3. Thanks to the assumed form of the power spectrum we may suppose that the 
realizations of the velocity field are P-a.s. analytic in the x variable, see e.g. I^. We will denote by 
E the mathematical expectation with respect to the measure P. 

The main results 

We consider the particle trajectory given as the solution to ()2.17() with the starting point x = and 
with the random field F(x) constructed in the previous section, and define the particle deviation 
from the mean position y{t) = X(t; 0) — vt. We also introduce 

Yi{t) := [ K[yi{s)fds, i = l,...,d. (3.1) 
Jo 

Suppose that the times Tg > are such that lim5^o+ Ts = +oo. Define the Cesaro limit of the mean 
square of the fiuctuation amplitude as 

C- hm E |y iTs)\^ := lim ^ Y (Ts) , (3.2) 

provided that the limit on the right hand side exists, whether it is finite, or not. 

We will distinguish two cases: as we have mentioned in the Introduction, when a G (1/2, 1) the 
Kubo- Taylor formula l|1.3|) diverges and we expect a behavior different from the usual diffusive limit. 
On the other hand, when a < 1/2 the diffusion coefficient given by H1.3() remains finite and the usual 
diffusive behavior would not be surprising. 
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The case a G (1/2, 1) 

We have argued formally in Section [21 that in this range the process y(t) converges on the time scale 
t ~ 0{5~'^^), H = l/(2a) to a fractional Brownian motion with the exponent a. Our first result 
confirms the predicted time-scale in this range of the parameter a. 

Theorem 3.1 Suppose that t > 0, a e (1/2,1) and p > 0. Let T/ := (^-2H(i+p)^i+p _ 
^-2H{i-p)^i-p ^ jjjidQj. ffig above assumptions about the random field F(x) we have 

C- lim E |y (r+) 1^ = +oo (3.3) 

and 



hm E|y(T5-)|' =0. (3.4) 



5^0 

The case a < 1/2 

In this case we expect that y(t) behaves diffusively on the time scale t ~ 0(6^'^), as in the situation 
when the correlation tensor decays rapidly in space, that is, when a is very negative. This time scale 
is confirmed by the next theorem. 

Theorem 3.2 Suppose that a < 1/2. Then, for arbitrary t,p > we have 

C- lim 6-PE |y (t^^) 1^ = +oo (3.5) 

and 

lim^E y (^tS-^^^-f'^^ =0. (3.6) 
Theorems 13. II and 13.21 are proved in Sections and El respectively. 



4 The proof of Theorem 13.11 
The lower bound 

We prove the lower bound (|3.3j) . This is done with the help of a general lower bound, which relates 
the long time behavior of the trajectory to the behavior of the resolvent near the border of the 
spectrum A = 0. In order to formulate it let us begin with some preliminary definitions. 

The basic spaces. For any p £ [l,-|-oo) and (j) ^ LP{^) we adopt the notation Dfc(/> := 
'^i4'{'Thek^)\h=o, where e^, k = 1, - ■ ■ ,dis the k-th vector of the canonical basis in M'^. The derivatives 
are understood in the sense. Let W^'"^ be the Banach space consisting of those elements (p E 
LP(F), for which 

ilH hid^'m 

We set 

C:= Pi WP'\ (4.1) 

l<p<+oo 

The spectral measure. The spectral theorem (see e.g. Theorem 1.4.2, p. 18 of j23j ) implies 
that there exists a complex vector valued spectral measure F(-) = {Fi {■),..., Fd{-)) defined over 
(M'^, ;S(M'^)) whose components take values in L^(P) such that 

F(x) = / e^''-'^F(dk). (4.2) 
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The spectral measure is Gaussian, tiiat is, for any set A € B{M. ) the 2d-dimensional random vector 
(ReF{A),lmF{A)) is Gaussian. Since the random field H4.2|) is real vector valued we must have 
F*(dk) = F{—dk). The corresponding structure measure equals 



K[Fi{dk)F*{dk')] = R^j{k)5{k - k')dkdk', i,j 



= l,...,d. 



(4.3) 



The Hermite polynomials. Suppose that "0 = (^i) ■ ■ ■ , 4^d), where G C^{M.'^), i = 1, . . . ,d 
are complex valued, even functions, that is, k) = ip*(k). We write 



Denote by H the subspace of (P) obtained by taking the closure of the linear span of the elements 
of the form (|4.4j) . It is a Gaussian Hilbert space in the sense of Definition 1.2 p. 4 of JB]- We can 
define then, see Definition 2.1 of ibid., the space of the n-th degree polynomials Vn as the L^-closure 
of all the elements of the form p{(j)i, . . . , (f)k), where p{-) is polynomial of at most n-th degree with real 
valued coefficients and (pi G H , i = 1, . . . ,k. It can easily be seen that Vn can be also characterized 
as the L^-closure of the space spanned by 



where ip : (R'^)" (C^)" is even, that is, 'ip{—ki, . . . , — k„) = ip*{ki, . . . ,k„), - see the Appendix 
for the definition of the multiple stochastic integral. For any n > we denote by if'"' := Vn Q Vn-i 
the space of the n~th degree Hermite polynomials, see Definition 2.1 of jl6j . Here V-i := {0}. It is 
well known, see Theorem 2.6 of ibid., that ^^(P) = 0^>q K''- 

A variational principle for the resolvent 

We define the random field V(x;a') := v + (5F(x;a;) and let X(t) be the solution of (|l.lj) with 
the initial condition x set at 0. The environment process is given by r][t) := T^{t-^ijo. It is an re- 
valued, deterministic, dynamical system, with P as its invariant measure. The corresponding group 
of Koopman operators = frjt{uj), t > extends then to a Co-continuous, unitary group on 

L^(P). Its generator L is given by 



Here V</> := {Di(f), . . . , Da(p). One can show, in the same way as it was done in see Lemma 4.1 of 
jl9j . that C is invariant under the action of the group (P*) and, since it is dense in L^(P), it is a core 
of the generator. 

Denote by Rxcp := (A — L)~^4) the resolvent operator defined for any A > and cp £ L'^(J?)- We 
can formulate now the variational principle that will be crucial for us in the sequel, cf. Lemma 2.1 




i=l 



(4.4) 




(4.5) 



L(P={v + 5F) ■ V<p, for (PeC. 



of 0. 

Proposition 4.1 For any f G 



(P) and X > we have 



(-Ra/,/)l2(P) = sup 2(/,(/>)i2(p) - -\\L(p\\l2(^p^ - A||</>||^2(p) ■■ (peC 



(4.6) 
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Proof. Let i?^ := (A — L)^ ^ be the symmetric part of the (bounded) operator Rx and let TZ be the 
range of Ri. As L is anti-selfadjoint, Ri is given exphcitly by 



1 



An elementary calculation shows that for any / G L'^(P) we have, as 7^ G D{L), 



{Rxf, f)LH¥) = {Rif, f)LH¥) = SUp[2(/, (/.)i2(p) - (0, (i?i) ^ 

1 



sup 



sup 



2(/,<A)l2( 
2(/, (/>)l2( 



A 



A 



((A-L)</., (A-L)0) 



\\\m 



2 

L2( 



>L2(p) : G 7^] 

G 7^ 



(4.7) 



(4.8) 



The last equality above uses anti-selfadjointness of L. As C is a core of L we can use it instead of IZ 
as the set of test functions in the variational principles ()4.8() - thus, (|4.6() follows. □ 

A lower bound on the variance of trajectory fluctuations 

The variational principle 1)4. 6p for the resolvent is used as follows. Let G = (Gi, . . . , Gd) G i>^(P), 
the next result concerns the lower bound of the second absolute moment of z(t) := G{7]{s))ds 
(here r]{t) is the environment process), cf. Lemma 2, p. 655 of ^H]. To formulate it we introduce 



^ ^ S / SI \ 

Z,{t) := ^ E[z^{s)fds = 2^ J I J {P'-^Gi,Gi)L^^)ds2 j dsi 



LO \0 



ds 



(4.9) 



and Z{t) = Ya=i Zi{t). Note that from we obtain, in particular that 

Zi{t) < i"'||Gj||^2(p-)/3. 



(4.10) 



The following proposition relates the small A behavior of the resolvent to the large time behavior of 
the process Z{t). 

Proposition 4.2 Let A > and define 



G{\) ■.= Y,{RxGi,Gi)L2^ 



(4.11) 



i=l 



Assume that t,P,p> are given. Let Ts := tS ^ . Then, there exist constants C^,5q > such that 

Z (r/+^) > aT|G(T,-i), V5 G (0,5o]. (4.12) 
Proof. Using 1)4. 9(1 and integration by parts we find that 



j e-^^Zi{t)dt = 2 jjjj e'^\P'^Gi,Gi)L2^^)dsids2dsdt 

0<S2<si<s<t 

+ 00 

= 2A-=' j e-^'^P''Gi,G^)L2(^^)ds2 = 2X-''{RxGi,Gi) 



(4.13) 
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On the other hand, since t Zi{t) is an increasing function, we can write that the utmost left hand 
side of (|4.13|1 is bounded from above as 



+00 A-(i+'') +00 



A-(i+p) 



(4.14) 



A-(1+P) 



A-(1+P) 



for some absolute constant C > 0. We have used H4.10() in the second inequality above. Performing 
now the integration on the utmost right hand side of (|4.14|) and recalling (|4.1I-{|) we obtain that 



r'a--^"''iir'.i|2 



3A4 



(4.15) 



Summing up over i we obtain (|4.12|) upon choosing A := Tj ^. □ 



The proof of the lower bound 

The lower bound in Theorem 13. II is proved as follows: note that 



y{t) = 5 fY{ri{s))ds 
Jo 



is of the form of the functionals considered in Proposition 14.21 fwith G = 5F). Therefore, to prove 
1)3. 3() we will first use suitable test functions in the variational principle 1)4. 6|) for the resolvent, and 
then use (|4.12j) to obtain the lower bound for Y{t) given by (|3.1|) . 

The test functions are chosen as follows: suppose that ^(k) G C^{M.'^) is a smooth complex 
valued, even vector function: ipil—k) = ■i/'?(k), i = l,...,d and denote by Vi the family of all 
random elements of the form 

(f){uj) = I V(k) • F{dk-uj). 



By virtue of (|4.6|) we have a lower bound 



{R\Fi,FilL-^ > sup 



2 

L2 



2 

L2 



Observe that 



2{Fi,4>)L2=2 



j ^(k) • r(k)f 



a{\k\)d\L 



|]^|2a+d-2 ■ 

where r(k) = [rij(k)]. Since we have 

Dj(i) = i j kjipik) ■ F{dk), 

the operator L acts as 



(4.16) 

(4.17) 
(4.18) 



(4.19) 



= (v + 6F) ■V(p = i 



V • k V(k) • F(dk) + 5 ip{k) ■ F{dk) k ■ F(dk') 



(4.20) 
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and 



we have 



2 

L2( 



I imm ■ m] (4-21) 



+5' 



2 _ /• (vk)2[r(k)V.(k)-V.(k)]a(|k|)dk 

LHP) - / |k|2a+d-2 



[r(k)V'(k) • ^(k)][r(k')k • k]a(|k|)a(|k'|)(ik(fk' 
(|k||k'|)2"+'^-2 

2 f f [ r(k)V;(k) •kq[r(kOk.^(kO]a(|k|)a(|k^|)rfkdk^ 

(|k||k'|)2"+'i-2 



(4.22) 



Using an elementary estimate 

(r(k)V'(k) • k') < (r(k)v^(k) • v(k))^/2(r(k)k' • k^/^ ^4 23) 

we conclude that 

{RxFi,Fi)L2 > sup[a*(0) : G Vi], (4.24) 
where 5j(0) is a quadratic functional given by 

W) := 2 1 [m . r(k)e.] - X I [mm . ^(k)] (4.25) 
1 /■ (v-k)2[r(k)V'(k) • V(k)]a(|k|)(ik 2(^2 r [r(k)V;(k) • V(k)][r(k')k • k]a(|k|)a(|k'|)dk(fk' 



Xj |k|2«+d-2 A 7 (|k||k'|)2"+rf-2 

Rd R2ci 

The maximizer of the functional given above equals 

V'*(k) =r(k)e,|A + ^ [(vk)2 + <52|k|27^]| , (4.26) 

r(k)ei-ei |j^,|2„+rf_2 -2|^l--j y ^j^TpTd^ > 0- 



where 



Thus, from (|4.24() we obtain 

d 

P{X) := Y,{RxFi,Fi)L2 > m{6), (4.27) 



i=l 

where 



miS) := (d - 1) / {a + i [(V. k)^ + 6^\M^n] }" ^|g£. (4.28) 

Rd 

We claim that m{5) > CA^~2a for < A < 1. To show this we will assume with no loss of 
generality that v = ei. The expression in 1)4. 28() is of the same order of magnitude as 

Adk ^ , /" (Ik , . s 



(A2 + A;2 + 2,52|k|2?t:)|k|2"+rf-2 - J (A2 + A;2+52|k|2)|k|2a+d-2 

lkl<l lkl<l 
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for some C > 0. Writing k = (/ci, 1) we can further transform the right hand side of (|4.29|) . It equals 

J [A2 + (1 + 52)A;2 + 52/2] (^2 + p)a+d/2-l ' ^ -^^l 

where ^ = |1|. Introducing the polar coordinates / = ^0036*, ki = Qsm.9 we can estimate (|4.3U() from 
below by 

} 7 cos-'^^ ededo )■ f dedo 

J J (A2 + £,2sin2 + ^2j2)^2a^l -^1 J J [A2 + ^2 (^2 + ^2)] ^2a-l " ^-^^ 




Substituting q' := X'^ gy/WTP we see that the right hand side of (|4.3ip equals 

7r/4 



de 







(^2+52) 



92+52)1/2a-i 

dg 







,2\ „2q!-1 



Since 2 — 2a < 1, the first integral above converges even for = 0, and the second also has a finite 



limit as A ^ 0. We conclude from the above and (|4.27j) that 

P(A) > CsA^-^"^ VA G (0, 1] (4.32) 

and some C2 > 0. 

With the lower bound (|4.32|) at hand we are ready to finish the proof of Theorem 13 .11 Let p' > 
be arbitrary and let be as in the statement of the theorem. We apply Proposition 14.21 with 
Ts := {T'^Y^^^^'^'\ Then, for sufficiently small 5q we have 

/ /\ lO^ 

for all 6 G (0, 60] and some C3, C4 > 0. Hence, 



— Y (r+) > C^d^Tg'^'''' = C452-2/f(l+p)(2a-p')/(l+p')^{l+P){2a-p')/(l+p') 

and 1)3. 3p follows, provided that p' > is chosen sufficiently small. 



The upper bound 

The upper bound 1)3. 4p in Theorem 13.11 is proved using an approximation of the trajectory by the 
correctors. 

The corrector fields 

Let fI^^'^'' := F and for any A > we define the corrector field of the first order in the direction of Cp 
Note that y^}] satisfies 



^P'^ ■= / A-ivk ^^^'^^^^ P = l,...,d. (4.33) 



(A - V • V)xj;l = p = l,...,d. 
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Let us define ^ = (F-J^, • • • , ^ij)^ where 



fSI ■= F ■ VxS = .// ^^i^ (4.34) 



p,X ■— ^ ^ Ap,A 

Then, we may write, using the definition of the first order corrector 



xp{t) -vpt = 6j Fpirj{s))ds = 5XJ xZ(vis))ds " ^ j v • VxZ(vis))ds. (4.35) 



On the other hand, we also have 

t 

4liv{t)) - xi;i(a;) = /(v + SF{r,{s))) ■ Vxj;i(r?(5))ds (4.36) 



t t 
= J^-yx^i{v{s))ds + 6 J F(2)(^(,)).Vxi;i(r/(s))ds. 



so that 1)4. 35(1 becomes 

t t 

xp{t) - V = I Fl,%{s))ds + 6\j X^^M^))ds + 5x^^M - ^yS^xm) (4.37) 



Now, we may iteratively define := F • Vx^xp and let x^xp be the solution of 

(A-v-V)xS = f(5. (4.38) 
Then, for any n > 1 we have a decomposition 

* n-l * n-1 

xp{t)-vpt = 5- / Fl%{s))ds + \Y,S' / xi'i(r?(^))d. + E^'[4!i(^)-4',A(^W)]- (4-39) 
'=1 '=1 

Technically, the most important results of this section are the following bounds for the correctors. 
Proposition 4.3 Suppose that a G (1/2, 1). Then, for each n > 1 we have 

||xl"^||L^<3^(l + |logA|)"/^ AG (0,1] (4.40) 

for some constant C > 0. When, on the other hand a < 1/2 we have 

llxf ||L^<3^(l + |logA|)"/^ AG (0,1]. (4.41) 

As a consequence we obtain the following. 
Corollary 4.4 For some constant C > we have 

||Fi")||L^<^(l + |logA|)"/^ AG (0,1] (4.42) 

for a G (1/2, 1) and 



\F^^^h^<^{l + \log\\r/\ AG(0,1] (4.43) 



for a < 1/2. 
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The proof of Proposition 14.^1 is rather technical and we postpone it until Sectional where the proof 
of Corollary 14.41 can also be found. 

In what follows we shall also need the following lemma which shows that the inner products of 
correctors with F are of a smaller order than one would naively expect from the Cauchy-Schwartz 
inequality and Proposition 14.31 

Lemma 4.5 Suppose that a < 1/2. Then, for a given n > 1 there exists a constant C > such 
that 



KxS,i^p)LH<3^(;^(l + |logA|) 
This lemma is also proved in Section 1^1 

The proof of the upper bound 

We now prove (|3.4|) . Using (|4.39j) we estimate 

n—1 c 



c 



(n+l)/2 



VA G (0,1]. 



(4.44) 



IE|y(Tr)r <C 



^^'"(T,-)^||Fi")||i.+j:<^^'(AT,-)^5;i|x«||i. + E^^T 



n-1 



1=1 



p=l 



1=1 p=l 



(4.45) 



Recall that Tj = (5~2//(i-p)^i p_ 'pj^jg together with estimates (fOU)) and imply that the right 

hand side of 1)4. 45(1 can be estimated by 



C 



n-1 



n-1 



1=1 



Choose A := 8'^^'^^^P\ we obtain then 

E |y (T^-) I' < Ct2(i-^) ( 5™P5-4H{i-p) ^ ^ 



21 



1=1 



n-1 



(4.46) 



1=1 



For an arbitrary /> > we can choose n sufficiently large so that the right hand side of 1)4. 46p is of 
order of magnitude o(l) so that 1)3.4(1 follows. □ 



5 The proof of Theorem 13.21 

For the most part the proof of the upper bound 1)3. 6|) in Theorem 13.21 is a repetition of what has 
been done in the corresponding situation in the previous case. Observe that H4.45|) still holds with 
T7 ■=t5-^^^-P\ We can choose now A := (T^ ) ^ and easily convince ourselves that ()3.6() holds. 

The heuristic reason why this direction is simpler once the upper bound in Theorem l3.1l has been 
obtained is that the time 6~'^ is much shorter that 5~'^^ , H = l/(2a) for a < 1/2. Therefore the 
fact that "nothing happens until the time 0((5~^)" is not very surprising in Theorem l3.21 The lower 
bound in Theorem 13.21 is more informative - it tells that something happens at the time 0(5~^) 
which is much earlier than the time scale 6~'^^ . 

We now prove the lower bound Let p' > and Ts := (t5-2)i/(i+p'). We shah further 

specify the parameter p' later on. By virtue of Proposition 14 . 21 we conclude that for a certain C > 
and sufficiently small 5 > 

y(tI+''')>C5^TIP{T,^), (5.1) 



16 



where -P(A) is given by H4.27() . Recall that the correctors x^\p ^-re the solutions of (|4.38|) for an 
arbitrary A > and N > n> 1. Let the remainders r[^J be given by 

[X-{. + 5F).V]r[l^=FZ-''\ (5.2) 

Multiplying both sides of (|5.2|) by r^^J and integrating out we get an obvious bound 

Observe also that 

TV 



n=l 

and therefore 



(iiAF„F,)^. > ix^F,),. -CY, n{xt:^'\F,)^A - CX-^l\5X-^l^ f{l + | logAI)^/^. 



n=l 



Lemma 14.51 allows us to obtain the following estimate 

N-l 

{RxFp, F,)l2 > {x% Fp)L2 - C ^ (5A-i/2)«| (1 + I log A|)"/2 - CX~^/\6X'^/Y {I + \ log A|)^/2. 

n=l 

(5.4) 

Choosing A := ^ we obtain from 1)5. ip and 1)5. 4() . with sufficiently large chosen, that 

j2-py (t5-2) > ct^/{i+p')s~P+^p'/(i+p')^ (5.5) 

which clearly implies (|3.5|1 . provided that — p + 4/9'/(l + p') < 0. □ 

6 The proofs of the corrector bounds 

In this section we prove the technical bounds on the correctors stated in ProDOsition l4.31 Corollarv l4.4l 
and Lemma 14.51 

The Feynman diagrams 

The proofs of the corrector bounds make an extensive use of the Feynman diagrams. Let us recall 
now the corresponding basic notions. Let Zn '■= {1, . . . yTi}, a Feynman diagram T (of order n > 
and rank r > 0) based on Z„ is a graph consisting of a set B{J^) C Zn of vertices from Zn, and a set 
E{J^) of e{T) edges connecting points in Z„ without common endpoints. So there are e{J^) pairs of 
vertices, each joined by an edge, and a{J-') := n — 2e{J-) unpaired vertices, called free vertices. An 
edge whose endpoints are m,n £ B is denoted by rrvh (unless otherwise specified, we always assume 
m < n). A diagram T is said to be based on B{T). Denote the set of free vertices by A{T), so 
A{T) = T \ E{T). The diagram is complete if A{T) is empty and incomplete, otherwise. Also for 
a given 1 < p < n let us denote by Ap{J^) the set of all free vertices that are less or equal to p 
and by Lp{J-) the union of Ap{T) and all the left vertices of edges mfi for which m < p < n. Let 
L{T) := Ln{T) denote the left vertices of all edges belonging to E{T). 
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Expressions for correctors 

We begin with some explicit expressions for the correctors Xm\- R-scall that the functions Xm A' 
m = 1, 2 . . . , d, are defined iteratively as the solutions of 

(A - V • V)xt\ = - F . Vx^:f, (6.1) 



with XmX^ solution of 



given explicitly by 

.(1) .. f 1 



(A - V • V)x^ A = ^m,A, m=l,...,d, 



^-,A ■■= I x-w-k = 1' • • • ' (6-2) 



Let us introduce auxiliary functions 



^n(ki, . . . , k„) 



for n > 1 and 

v-vn— 1 1 

A - iv-ki A - iv-(2^;^iki) 

Then, a simple induction argument shows that the n-th order corrector in the direction of the vector 
Bp is given by 

■= / h^r\^i) ® /^2(ki,k2) ... hniki, . . . ,k„) •F(dki) F(dk2) ... F(dk„) (6.3) 

and the fields fJ"^ defined in dHU) are fJ"^ = (F^^"\ • • • , ^i?), where 

:=i""2ye„®/ii(ki) 0...® /i„_i(ki,...,k„_i)-F(dki) ® F(dk2) . . . F(dk„). (6.4) 

The proof of Proposition 14.31 
The Feynman diagrams expansion 

We now prove Proposition 14.31 an L^-bound for the correctors. Using the Feynman diagram expan- 
sion (|A.7|) of a multiple stochastic integral from the Appendix we can write that 



X?=Y.^^x\n (6.5) 



where the summation extends over all Feynman diagrams J- of the set Z„. For a given diagram 
with free vertices A{J^) = {ni, . . . , Ua} and edges E{J^) = {pq} we set 

xI:?a(.^) /• • •/ r'\K„. . . ,k„J • F(dk„J ^...^ FidKJ :, 
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where : • : denotes the orthogonal projection onto the space of a-th degree Hermite polynomials 
H-"-, cf Definition 3.1 of JL^, and f""'^ : {M.'^)'' {R'^)" is given by 



Using Theorem 3.9 p. 26 of one concludes that 



„ J 1 



a(|k„J) 



|2a+d-2 



(6.6) 



where the summation extends over all permutations vr : A{T^ ^(•^)- We assume with no loss of 
generality that v = ei. Also we write kj = (/cj, Ij) G R x M'^^-'^, splitting out the first component of 
the vector k, . Then, with this notation we have 



K K 



ic....„.(k„,,...,k„ji< n — ^ 



K -K 



A - i E 



n 



pqeE{r) 



r<"nl 



X < 



L |ipl<^ 



yr 6{lp + lq)dlpdlg 



The last integral can be estimated by CY[pq^E{j^)(^ + l^pl) conclude that the right hand 

side of (|6.6p can be estimated by 



.(^)!/ir-Y(k„,,...,k„jn*^ 



p=l I "Pi 



(6.7) 



K K 



< 



-K -K 



where 



K K 



n 

n — - 



-K -K 



A - i E ^« 



1=1 



n 



pq&E{T) 



[l+\0g+\kp\)5{kp+kq) 
\h |2a-l 



dkpdkq . 



The proof of ^^^^ 

Suppose first that a G (1/2,1) - we need to prove the estimate ()4.40() from Proposition 14.31 Note 
that after changing variables kj = Xkj and dropping the primes we get the estimate of the left hand 
side of (|6.7|) as 

\\xt\i^)\\h<CX-'^-{l + \logX\ri, 
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where 
with 



gikrix, ■ ■ ■ ,kna) :— f ... f Y\ ~ JJ I J pa-l ^ —dkpdkq. 

-oo -oo i=l 1 + I ^ A;,| P9G£{^) ' 

/=1 

We shah show by an elementary calculation that / < +00 - this is all that remains to prove the 
estimate (|4.40p . We present the details for the convenience of the reader: note that 



!=[...[ fr — X ft L — X fr 

i=l \ 2^1=1 '^H j'=l ^ ^ \ 2^l'=l'^l'\ m=l 



\ (l + log+|fe„„|)dfe„„ 



X 



k |2a-l 

{I +lOg+ \kp\)6{kp + kg)dkpdkg yr (1 + 10g+ \kp,\)S{kp, + kg,)dk'j^,dkg 



[l+\0g+\knj)dkn„ 



kry 



\2a-l 



/lb I L ^ Li 

TT — X TT — X TT 
^„ j=l \ 2^1=1 j'=l + I Z^/'=l I m=l 

-pr (l + log+|A:p|)dfcp -py (l + log+IVDdfc;, 

^ li It |2a-l X J_J_ |, |2q:-1 

We adopt the rule that in the sums above k^ + kg = and /c(, + /cl, = 0, if fs is an edge. Now, the 
integration over \k\ < 1 is not a problem so in order to show that / < +00 it suffices only to prove 
that 

n ^ n ^ 

i'= TT — xTT — 

^ W (l + log+|A:„„|)cifc„^ (f + bg+M^^ TT (1 + log+lVlK, . , 

^ ii l + |2"-i ii l + |A:J2"-i ii l + 

Using Holder inequality with p such that p(2q — 1) = 2 and = 2/(3— 2a) > 1 such that = 1 

leads to /' < C(//)V9, with 



/ n ^ ^1 

IT i+iy-p n X n X J] X J] dK,. 

R„ p=l I ^^=1 'I p'=l I Z^/'=l '^'F I rGAm mgim r'eLm 



Now, we introduce new variables rji, . . . ,r]n as follows: first, we look at the sums 5m = X^jLi % with 
m = 1, . . . , n, with the terms A;,. + fe^ = if fs is an edge. We say that 771 = fei, then we take j2 
that is the smallest number larger than one which is either free or a left vertex, and let 7/2 = Sj^. 
Note that either r]2 = k2 if k2 is not a right end, or 772 = ^3 if (^1^2) is an edge. We continue in 
the same way: having defined r]i_i = Sj^_^ we take ji to be the first index larger than which 
is not a right end of an edge, and take rji = Sj^. In this way we will pick n — e{T) = a{J^) + e{J^) 
sums out of 5m, m = 1, . . . , n and will define rii,r]2, ■ ■ ■ , 'i]n-e{:F)- Note that the change of variables is 
lower-triangular. Then, we do a similar procedure with the /c'-variables - we look at Pm = SJLi 
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and set rjn-e+i = Pqi, I = 1, • • • jC, where gi < . . . < ge are the left vertices. This defines ah rji, 
I = l,...,n and keeps the change of variables lower-triangular. Hence, it is invertible with the 
Jacobian equal to 1. As a consequence, we may drop all sums in the denominator which are not 
equal to one of r]j and we have 

. n ^ 

It follows that / < +00 and the proof of (|4.4nj) is complete. 



The proof of (li:iT|) 

It remains to prove the estimate 1)4. 41() in Proposition 14.31 Suppose now that a < 1/2. The right 
hand side of (|6.6j) can be estimated by 

K K ^ 
I '■= C J ... J {kn-^ , ■ ■ ■ , kua ) dkrip , 



(6. 



where 



-K -K 
K K 



p=l 



f\{kni 1 ■ 



We have 



KC 



k 

1 '^ria 



K 



n 



n 



-K -K 



K n 



1=1 



6{kp + kq)dkpdkq. 



pq(iE{T) 



n 



n 



(iA;„^ X 5{kp + kq)dkpdkq X ]^ S{kp/ + kg/)dk'p,dkg 



m=l 



We introduce the new variables 771, . . . , r/„ as in the previous case. With the help of this transforma- 
tion and replacing the remaining n of the denominators above by 1/A we conclude that 



^ J-nK J~nK J Jl^+ Wjl ^ 

for some C > and ah A G (0, 1]. Thus (|OT|) follows. □ 



rnK n 



log ( i ) + 1 



(6.9) 



The proof of Corollary 14.41 

We now prove the estimate 1)4. 43() in Corollarv 14.41 Using the Cauchy- Schwartz inequality we obtain 

l|Fi")||i.<||F|U4||Vx?k- 
However, since all the norms, 1 < p < +cxd on Vn are equivalent, see Theorem 5.10 p. 62 of jl6j . 

(n) 

and F G 7^1, G we conclude that 

l|F$:)|lL^<C||F||^.||Vxf 11^. 

for some constant C > 0. Thanks to the fact that the spectral measure F has a compact support 
the gradient operator V restricted to Vn is bounded (this can be seen by a direct calculation Dp on 
elements of the form (|4.5|) ). This in turn implies (|4.43jl . □ 
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The proof of Lemma 14.51 

Finally, we prove the estimate (|4.44)) in Lemma 14.51 We only need to be concerned with n odd, for 
otherwise the left hand side of H4.44() vanishes. Using expansion ()6.5() we can write that 



(xS'^p)i^ = T.(^x],i^)^^p)L- (6.10) 



and the summation extends over those Feynman diagrams for which a{J^) = 1. Suppose that q is 
the only free vertex of J^. Then, 



X |],"|2a+i-2 ^(^1 + kn+l)rfk,dk„+i J] Ti,,. (k.) "j^^^^ (5(k. + k,)dk,dk, . 

Using the change of variables made in the course of the proof of Proposition l4.3l we conclude therefore 
that 

Kxg(-^)>^p)LH<^(i + iiogAir+\ 

which yields ()4.44p because 2e + 1 = n. □ 

A Multiple stochastic integration 

Suppose that P is a Gaussian, homogeneous, Borel measure over the space il, introduced in Section 
131 Denote by F the corresponding Gaussian vector valued spectral measure on (M*^, ;S(M'^)). We 
suppose that the structure measure R((ik) has the density r(k) = [rjj(k)] w.r.t. the Lebesgue 
measure. Let r(k) :=tr r(k) V 1. For a given integer n > 1 we consider the Borel measure on (M'^)^"' 

2n 

M2„(dki, . . . , dk2n) = H ^^^P + n ^(ki)dkj-. 

y^eSn pqeT i=i 

By £^ we denote the completion of the space of all complex valued, bounded Borel measurable 
functions ip : (M'^)" — > (C^)" in the norm 

£2 := j ■■■ j (|V'(ki,...,k„)p + |^(k„+i,...,k2„)P)M2„(dki,...,dk2„). 



Suppose that i G Z'^ and > 1. We define 

□7v(0 := [k e M'^ : < kj < 2'^ {ij + 1), Vj = 1, . . . ,d]. (A.l) 

For i := (ii, . . . , i„) G [Z'^Y we let 

□jv(i) :=njv(n) X ... X njvlin). (A.2) 

By T) we denote the family of all such boxes. Its subfamily n is called an admissible dyadic partition 
of (IR'^)2" if 

(PI) Ungnn = (]R')'"> 
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(P2) for any two boxes □ 7^ □' € 11 we have □ n = 0, 
(P3) there exists do > such that |n| > do for all Dsn. 

A set function c : 11 — > (C'')"' is called admissible if c(n) = for all but finitely many boxes from 
an admissible 11. We denote by A the family of all admissible set functions. 

For any multi-index j := (ji, . . . ,jn) £ {1, . . . , d}" and □'^^^(i) given by (|A.2|) we let 

n 

Fjp^(i)] = n4p^(v)]. (A.3) 

p=i 

Suppose that c : II (C^)" is admissible. We define then V'(k) := c(n) for all k G □. With 
some abuse of notation we call such a function admissible and denote by An C £^ the space of all 
admissible functions. For any admissible function ^ E An we define the n-tuple stochastic integral 
letting 

□en j 

We shall also write /. . ./ ■ ■ ■ , k„)F((iki) ... F((ik„) to denote Below, we list some 

of the properties of T(^). They are elementary and their verification relies on the application of the 
definition so we leave this task to a reader. 

Proposition A.l (i) An is dense in £^ in the norm \\ ■ \\c'^- 

(a) The stochastic integral given by (|A.4I) is well defined, i.e. if there exist two admissible set 
functions c\ , C2 corresponding to a given ip then the respective definitions of the stochastic integrals 
are identical. 

(Hi) We have T(ai^/'i + a2il^2) = + o,2'^{ip2) 

(iv) Suppose that V'l, . . . , V'n £ Ai. Then V'l (X" • • • (X" £ An and 

n 

(v) We have 

E[T(^«)r(^(2))] = X; /•••/ ^£...Jki,...,k„)(^(;)^^,...,,,J*(k.+i,. 



2n 



X H Rpg{kp)6{kp - kg) II dkj. (A.5) 

pqGJP j = l 

As a direct consequence of property (v) and the definition of the norm on £„ we obtain. 

Corollary A. 2 An 3 ip 1-^ ^-^ cl continuous functional that extends to the entire The 

extension shall be called an n-tuple stochastic integral w.r.t. the spectral measure F{dk). It has 
properties (iii)-(v) from Proposition \A . 1\ 

Proposition A.3 For any Tp £ C'^ we let 

U'^^p{\^l,... ,'kn) :=exp|i ^^kpj • x| ^(ki, . . . , k„). 

Then, for any x € we have U^xp € £^ and 

I(i;){T^uj) =I{U''ip){uj), F-a.s. (A.6) 
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Proof. We have 

for all x,y e R'^. Hence F((ik;Tya;) = e'^ yF{d'k;uj) and (fX!6|l follows. □ 

Remark A different construction of a multiple stochastic integral can be obtained using the 
approach of Chapter 7.2 of UH]- Denote the stochastic integral, defined there, by 



J. . .J V(ki, . . . , k„) • F(dki) ... F(cik„). 



Let J'-" be a Feynman diagram labelled by {!,... , n}. Let ni < . . . < Ua he all the free vertices of 
and let pq £ E{T) denote the remaining e := (n - a)/2 edges. Set T{T)'il) : [W^f {C^Y 



) := •••y V'ii...i„(ki,...,k„) ripi^(kp)5(kp + kg)(ikpdkg 

ipiq pq&E{T) 

for ailing , . . . , in„ = 1, . . . , d. Let 

:= y. . .j^'^^ T{T)^{kn,,. . . ,k„jF(dk„J ... F(dk„J. 
Using Theorem 7. 25 p. 99 and Corollary 3. 17 p. 28 of ibid, one can conclude that 

j ...j ^(ki,...,k„)F(dki)®...®F(dk„) = ^/(.F), (A.7) 

where the summation extends over all Feynman diagrams labelled by {1, . . . , n}. We could use then 
equality ()A.7|) as the definition of the multiple stochastic integral. 
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